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One-electron self-interaction and an incorrect asymptotic behavior of the Kohn-Sham exchange- 
correlation potential are among the most prominent limitations of many present-day density func¬ 
tionals. However, a one-electron self-interaction-free energy does not necessarily lead to the correct 
long-range potential. This is here shown explicitly for local hybrid functionals. Furthermore, care¬ 
fully studying the ratio of the von Weizsacker kinetic energy density to the (positive) Kohn-Sham 
kinetic energy density, rw/r, reveals that this ratio, which frequently serves as an iso-orbital in¬ 
dicator and is used to eliminate one-electron self-interaction effects in meta-generalized-gradient 
approximations and local hybrid functionals, can fail to approach its expected value in the vicin¬ 
ity of orbital nodal planes. This perspective article suggests that the nature and consequences of 
one-electron self-interaction and some of the strategies for its correction need to be reconsidered. 


I. DENSITY FUNCTIONAL 
APPROXIMATIONS AND THEIR KOHN-SHAM 
POTENTIALS 

During the past decades, Kohn-Sham density- 
functional theory (DFT) [TJ [5] evolved into a stan¬ 
dard tool for electronic structure calculations of atoms, 
molecules and solids. The decisive quantity of DFT is the 
exchange-correlation (xc) energy functional, E^c, which 
contains all electronic interaction beyond the classical 
electrostatic Hartree contribution, E^. ifxc in practice 
has to be approximated, and the approximation used 
governs the accuracy of a DFT calculation [1U. It 
is one of the puzzles of DFT that explicit density func¬ 
tionals such as the generalized gradient approximations 
(GGAs) can predict binding energies and bond lengths of 
complex many-electron systems reliably, but make sub¬ 
stantial errors in describing simple one-electron systems. 
The underlying problem is well known as the one-electron 
“self-interaction problem” [^r For the exact functional, 
Exc + Eh will vanish for any one-electron ground-state 
density because one electron does not interact with itself 
- but most approximate functionals yield a spurious fi¬ 
nite value for this case. Following Ref. [S] a functional 
is considered to be one-electron self-interaction free if it 
fulfills the condition 

Enlnia] + Exc[nia] = 0 , ( 1 ) 

where rua = |(^i(T(r)p designates a single spin-orbital 
density. 

Self-interaction plays a decisive (although not the only) 
role in the (un)reliability of density functional theory 
calculations, and its consequences are particularly pro¬ 
nounced, e.g., in questions of orbital localization ISHS], 
ionization processes [Ml], charge transfer HMSj, and 
for the interpretability of eigenvalues and orbitals, e.g., 
as photoemission observables [T6H23|. 

Many of these observables can also be directly related 
to properties of the Kohn-Sham exchange-correlation po¬ 
tential, which is defined as the functional derivative of the 


xc energy with respect to the ground-state density n(r), 
i.e., Uxc(i’) = It is generally expected that there is 

a close relation between freedom from self-interaction and 
xc potential features. The field-counteracting term that 
is important for obtaining correct response properties is 
one such feature [13] US] . Another example, and proba¬ 
bly the most prominent one, is the long-range asymptotic 
behavior of the xc potential [261 [27] , 


(Hartree units are used here and throughout.) In this 
perspective article we focus exclusively on Kohn-Sham 
theory, i.e., on a local multiplicative xc potential, as op¬ 
posed to orbital-specific (non-multiplicative) potentials 
that arise in generalized Kohn-Sham theory |28j and are 
used in the standard application of hybrid functionals. 
In the Kohn-Sham approach, the local xc potential mod¬ 
els the interaction of one particle with all others and it 
therefore appears intuitively plausible that a functional 
that is not self-interaction-free cannot show the correct 

— 1/r long-range asymptotic behavior: As one particle 
of a finite, overall electrically neutral systems ventures 
out to infinity, it will “feel” the hole of charge 1 that it 
left behind in the total charge. This gives rise to the 

— 1/r potential asymptotics (see, e.g.. Ref. |3|, p. 242 for 
a more detailed argument along these lines). However, 
a particle that spuriously self-interacts will “feel” itself, 
and thus not the proper hole. Consequently, the poten¬ 
tial will not have the proper long-range decay. 

The correct asymptotics of the xc potential has proven 
to be important for a variety of physical quantities. It 
plays a prominent role for obtaining stable anions in 
DFT, it leads to a Rydberg series in the Kohn-Sham 
eigenvalues and generally to unoccupied eigenvalues of 
improved interpretability, and as a consequence allows 
for improved accuracy in the prediction of various re¬ 
sponse properties [29H33|. The correct asymptotic be¬ 
havior is also important for the ionization potential (IP) 
theorem milTlIMlISSl, which states that the negative of 
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the highest occupied Kohn-Sham eigenvalue —eho should 
correspond to the vertical IP, and for developing func¬ 
tionals that allow for approximately predicting IPs from 
ground-state eigenvalues I3S1EZ1- 

There have been fruitful attempts to incorporate the 
correct behavior in the limit |r| —^ oo directly into the xc 
potential |38H41j , leading to improvements in the descrip¬ 
tion of some of the aforementioned properties. However, 
since directly designed potential expressions are typically 
not functional derivatives of any energy functional, the 
use of such “potential only” approximations is necessarily 
limited, as discussed, e.g., in detail in Refs. [HHH]. 

A functional that combines freedom from self¬ 
interaction and the correct asymptotics of the potential 
is exact exchange (EXX), being defined as the Fock inte¬ 
gral evaluated using Kohn-Sham orbitals tpia (r), where i 
labels orbitals in spin channel u: 

= E 

i,j — l 

(3) 

Here, is the number of electrons with spin cr. Treat¬ 
ing exchange exactly with a local Kohn-Sham potential 
leads to a significant improvement in the quality of Kohn- 
Sham eigenvalues when comparing to (semi-)local func¬ 
tionals [ini 113 SS] ■ EXX also tends to increase Kohn- 
Sham gaps [3 ITTHSO] , leads to a desired particle number 
discontinuity in static [51] and time-dependent [52] sit¬ 
uations, and improves the description of charge transfer 
[23 [2S] , dissociation [^ and ionization processes [S2] . 

However, using bare EXX is known for its rather poor 
description of binding energies and structural properties 
(see, e.g.. Refs. [n3[5S], and Ref. [3, chapter 2). Adding a 
(semi-)local correlation term to EXX hardly improves the 
situation and typically leads to results that are inferior 
to the ones from (semi-)local functionals. The reason 
for this failure is the well-known incompatibility of the 
fully non-local Fock exchange with a purely (semi-)local 
correlation term |5hj . 

A class of approximations which has been designed to 
remedy this incompatibility is the one of local hybrid 
functionals [53 EH], sometimes also called hyper-GGAs 
[56] ■ Whereas global hybrid functionals [591 - 163 ^ 

constant, fixed fraction of Fock exchange with (semi- 
)local exchange and correlation, local hybrids replace the 
fixed fraction by a density dependent local mixing func¬ 
tion (LMF). Both types of hybrids originate from the 
concept of the coupling-constant integration, i.e., adi¬ 
abatic connection scheme [55] EH]. Global hybrids are 
successful in modeling the coupling-constant averaged, 
integrated energy. Local hybrids can go one step fur¬ 
ther and aim to model the coupling-constant curve itself 
[55] instead of just the integral. Thus, in contrast to the 
global hybrid functionals which are used in practical ap¬ 
plications of DFT and combine GGA components with 
about 25 % of exact exchange, local hybrids can incor¬ 
porate full exact exchange and can be fully one-electron 


self-interaction-free. 

An early local hybrid with reduced one-electron self¬ 
interaction error showed promising results for dissocia¬ 
tion curves and reaction barriers, but its accuracy for 
binding energies was limited [58]. A self-consistent im¬ 
plementation of a local hybrid functional was given in 
Ref. [ 53 , and over the years several local hybrids were 
constructed, using different LMFs and (semi-)local ex¬ 
change and correlation functionals [55H75] . striving to 
reach greater accuracy by refining the position-dependent 
mixing of nonlocal and local components. Many of these 
functionals rely on the concept of an iso-orbital indicator, 
i.e., a functional that allows one to distinguish regions of 
space in which the density is dominated by one orbital 
shape from regions of space where several orbitals of dif¬ 
ferent shape contribute to the density. The most promi¬ 
nent iso-orbital indicator, which goes back to a long tradi¬ 
tion of using kinetic energy densities in density functional 
construction [Z3[ZS], is the ratio of the von Weizsacker 
kinetic energy density rw to the positive (as opposed to 
other possible definitions, see, e.g.. Ref. [55]) Kohn-Sham 
kinetic energy density r, discussed in detail below. 

By using full EXX and an iso-orbital indicator, lo¬ 
cal hybrids aim at being one electron self-interaction-free 
and producing a Kohn-Sham potential with the proper 
long-range asymptotic decay. They are a paradigm class 
of functionals designed for simultaneously curing both 
of these two prominent problems of (semi-)local density 
functionals. In the following, we therefore use the ex¬ 
ample of a local hybrid functional to shed light on the 
relation between a functional’s self-interaction and its 
potential asymptotics, as well as the properties of the 
tw/t indicator. We argue that quite generally a one- 
electron self-interaction-free energy does not guarantee 
the correct long-range potential, and that tw/t loses its 
indicator ability in the vicinity of nodal planes of the 
highest-occupied molecular orbital (HOMO). 

II. CORRELATION COMPATIBLE WITH 

EXACT EXCHANGE: THE LOCAL HYBRID 
APPROACH 

The xc energy functional can be written as 

Exc[n] = j n(r)exc([«];i')d^r, (4) 

with exc([fT-];r) denoting the xc energy density per parti¬ 
cle. The definition of exc(r) is not unique and subject to 
a gauge-dependence [72] ■ Yet, for local hybrid function¬ 
als it has become common to define this energy in the 
form 

e))L(r) = er(r) + /(r)(e^'(r) - e^ir)) + efir). (5) 

Here, e®’’'(r) marks the exchange energy density per par¬ 
ticle corresponding to the EXX energy of Eq. (§. This 
nonlocal term is mixed with (semi-)local exchange and 
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correlation energy densities e®'(r) and e®'(r), respectively. 
The position dependent mixing ratio /(r), which is itself 
a density functional, marks the LMF. 

Often, the LMF is designed in a way that aims at elim¬ 
inating the one-electron self-interaction error of Eq. Q 
that is inherent in most density functionals. An estab¬ 
lished method for reducing self-interaction effects is to 
detect regions of space where a single Kohn-Sham orbital 
shape dominates the density (“iso-orbital regions”), and 
then enforce E q. ([l|) in these regions. One of the most 
popular [5H1 Ezi Esnzii m iza iHM indicator functions 
for detecting iso-orbital regions is 


ffW = 


r(r) 


( 6 ) 


where rw(r) = |Vn(r)p/(8n(r)) denotes the von 

Weizsacker kinetic energy density and r(r) = 
i Sfci is the positive Kohn-Sham kinetic 

energy density. In iso-orbital regions, r(r) —)■ rvi/(r) and 
therefore g(r) —>■ 1. In the case of a slowly varying den¬ 
sity, T^(r) —0 and, since T(r) remains finite, g{r) —>• 0. 
This indicator function is typically a decisive ingredient 
in the LMF, /(r), of local hybrids. With its help one 
can construct /(r) such that Eq. ([^ reduces to correct 
limiting cases, e.g., e®'(r)-|-e®'(r) for slowly varying densi¬ 
ties, and e®^(r) for single orbital regions. The latter case 
additionally requires that e®^(r) vanishes in single-orbital 
regions, a condition that we discuss below. 

In the asymptotic limit, |r| —?> oo, the xc energy density 
for a finite system should be dominated by e“(r). When 
e®^(r) vanishes sufficiently fast in the asymptotic region 
(a condition that is usually fulfilled), then 


lim /(r) = 0 (7) 

|r|—>-oo 


is the requirement that one aims at, because it leads to 
the correct asymptotic limit of the xc energy density per 
particle 

( 8 ) 

|r|-).oo 2r 


(Note the difference to the asymptotic limit of the xc 
potential, see Ref. [31]). 

Since for a finite system each Kohn-Sham orbital de¬ 
cays exponentially with an exponent set by its eigen¬ 
value [M], the density is asymptotically dominated by 
the HOMO density, i.e., becomes of iso-orbital charac¬ 
ter. Therefore, g{r) can be used in the construction of 
the LMF to realize Eq. Q. 

Considerations of the type discussed above are inherent 
to many density functional constructions. As a particu¬ 
lar example for a local hybrid functional we here use a 
recently proposed, physically motivated LMF [55] , which 
reads 


ft{r) 


l-^C^(r) 

1 + c<2(r) 


(9) 


The function g{r) in the numerator is multiplied by the 
squared spin polarization C(r) = (n^(r) —n^(r))/(n^(r)-|- 
n^(r)), which lets the LMF not only identify iso-orbital 
regions, but also correctly distinguish between true one- 
orbital regions, and regions with two identical spin- 
orbitals. The function g{r) is used in such a way that 
/t(r) vanishes for one-orbital regions, as required. The 
use of the reduced density gradient 


^ ^ Vs/ 16$2(((r))’ n7/3(r) ’ 


( 10 ) 


where oq is the Bohr radius and $(C(r)) = 

in the denominator of /t(r), 
ensures the correct behavior of Exc under uniform coor¬ 
dinate scaling r —>• yr [SSJ [53]. The density transforms as 
n^{v) = and as a consequence Eq. (j^ uses full 

exact exchange in the sense of 


Ex 


lim 

7->'Oo E^[n^] 


= 1 . 


( 11 ) 


The function t^{v) is multiplied by a parameter c that we 
cannot determine, at least presently, from fundamental 
constraints. It allows for adjustments in the functional 
ansatz. In the case of slowly varying densities, /t(r) —>■ 1 
and Eq. ([^ reduces to its purely (semi-)local compo¬ 
nents. As an aside we note that this LMF comprises the 
one of Ref. [SS] as the special case c = 0 and C(r) = 1 Vr. 
We denote this case by fo{r), i.e., /o(i’) = 1 — Y(r)^ ■ 

For the semi-local exchange we use the 

LSDA [57], i.e., ex(r) = ex®°^(r), whereas 

e®^(r) = e);'®°"^(r). The additional 

multiplication with the numerator of Eq. consistently 
reduces Eq. ® to pure EXX in the one-spin-orbital 
case, where e^^^(r) alone does not vanish. 

The general questions that we discuss in this perspec¬ 
tive article, i.e., whether there is a relation between self¬ 
interaction and the xc potential asymptotics and in how 
far the iso-orbital indicator tw /t can be used to enforce 
freedom from self-interaction, can be scrutinzed with the 
local hybrid of Eq. © as an instructive example. 


III. THE KOHN-SHAM 
EXCHANGE-CORRELATION POTENTIAL OF 
LOCAL HYBRID FUNCTIONALS 

In order to implement local hybrids self-consistently 
within the Kohn-Sham scheme, one has to find the lo¬ 
cal multiplicative xc potential corresponding to the en¬ 
ergy of Eqs. Q and ([^. The fact that local hybrids 
use EXX and typically also T(r) makes them explicitly 
orbital-dependent. Therefore, the local xc potential must 
be obtained from the optimized effective potential (OEP) 
equation (see, e.g., [45l [55l l88ll9T] L The computational 
effort can be reduced significantly by employing the ap¬ 
proximation of Krieger, Li and lafrate (KLI) [^. For the 
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local hybrid of Eq. it has been shown that the total 
energy Etot and the highest occupied Kohn-Sham eigen¬ 
value Eho obtained with the KLI approximation agree 
quite well with the ones from the full OEP [S3]. Fur¬ 
thermore, it is a general finding[53] that the KLI ap¬ 
proximation does not affect the potential asymptotics to 
leading order. In the actual calculations presented in the 
following we therefore always use the KLI approximation. 

In the OEP (and KLI) scheme the chain rule for func¬ 
tional derivatives [IS] relates the derivative with respect 
to the density to the derivatives with respect to the or¬ 
bitals, 




1 (5Exc[{yj}] 


( 12 ) 


By defining P(r) = n(r) (r) and (5(r) = n(r) e®’'(r) 
one can write 


= 


Sip^air) 


J f{v') (P(r') - Q(r')) d^r' 




( Sfjr^) 


P(r') d' 




' f Sf{r') 


\6ipia{r) 


Q{r')d\' 


-fm 


( SQjr') 

U‘/5*<T(r) 


A^r' 


(17) 


From the structure of the OEP equation it further follows 
that to hrst order 


lim Uxca(r) = , lim ujv„o-(r), (13) 

|r|—>-oo |r|—>-oo 


i.e., the functional derivative with respect to the HOMO 
in general determines the potential asymptotics |84j . 
Therefore, investigating the HOMO functional deriva¬ 
tive is the key to determining the asymptotic behavior 
of an orbital dependent functional’s xc potential. When 
one takes the functional derivative (with respect to the 
orbital) of a local hybrid one obtains three terms, corre¬ 
sponding to the three addends in Eq. (§: 

ut{v) = <r(r) + + ulT\r) (14) 

Evaluating the asymptotical behavior of each of these 
three terms for the highest occupied orbital allows one 
to predict the potential asymptotics. 

The first term can be derived directly from Eq. ^ and 
reads 







(15) 

This term evaluated for the HOMO indeed provides the 
correct asymptotic behavior [J^ 


The first two terms consist of (semi-)local components 
and thus vanish exponentially. Evaluating the third term 
on the other hand is not as trivial as it contains the non¬ 
local quantity Q(r) as well the functional derivative of 
the LMF with respect to the corresponding Kohn-Sham 
orbital. For the LMFs addressed in this perspective we 
find that this term does not contribute to the asymp¬ 
totics either (see Ref. [83] and Appendix |B for details). 
Thus, only the fourth term in Eq. ( [17] ) is relevant in the 
asymptotic limit. 

Thus, only the fourth term in Eq. 0 is relevant in 
the asymptotic limit and therefore 


c—nl 


(r) 


2‘^L(r) 


N, 








i=i 


r — r 


i=i 




r - r' 


■ (18) 


The first term in this equation equals of 

Eq. ( [I^ , locally multiplied by /(r)/2. Due to Eq. Q 
it vanishes faster than the leading term of v}^, which is 
given in Eq. (16). 

The second term, however, is of a different structure, as 
it evaluates the LMF under the integral. By considering 
the HOMO level, its asymptotic limit is 


u 


exx 

N„<7 


(r) 




(16) 


c—nl 


.(r) 


|r|—>-oo 






r - r' 


d^r'. (19) 


The third term ^^^.’^'(r), on the other hand, does not 


contribute to the asymptotics of Eq. (16) as it decays 


exponentially due to its purely (semi-)local nature. 

Evaluating the second term on the right-hand side of 
Eq. (14) requires careful consideration. Intuitively, one 


might expect that an asymptotically vanishing LMF will 
surpress any asymptotic contribution of this term to the 
potential. In the following we check this expectation. 
Details of the underlying calculation for both LMFs used 
in this work, i.e., /t(r) and /o(i’), can be found in Ref. [53] 
and in Appendix [^ Eq. (Bl), respectively. 


This corresponds to a Hartree-like potential caused by 
the spin-orbital density of the HOMO averaged over all 
space, with the LMF as a weighting function. Thus, this 
term gives a Hnite contribution in the asymptotic limit 
despite of Eq. ([^. 

Now, when adding the asymptotically significant com¬ 
ponents, Eq. 0 and Eq. ( [I^ , for the evaluation of 
Eq. (13), we arrive at 


(r) 


|r|—>-oo 


( 20 ) 
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FIG. 1: The xc potential i>xct(r) for the C atom along the 
x-axis, (see Appendix [A| for definition), computed using ft(j) 
with c = 0. 5. Also displayed is the asymptotic curve according 
to Eq. (201 with 7 t(c = 0.5) = 0.716 and the correct asymp¬ 
totic — 1/r. The inset shows the potential plotted along the 
complete axis. 


FIG. 2: xc energy density per particle exc(r) for the C atom 
along the x-axis (see Appendix for definition), computed 
using /t(r) with c = 0.5. Also displayed is the corresponding 
asymptotic slope of — l/(2r). The inset shows the energy 
density plotted along the complete axis. 


Here, the parameter 7 ^ denotes the reduced slope of 
the potential asymptotics, which can numerically be ex¬ 
tracted from a self-consistent Kohn-Sham calculcation 


via 


7<t = 1 - ^ J dV. (21) 


Eq. (21) is a central result of this work, as it demon¬ 
strates that a local hybrid of the form of Eq. ([^ does not 
lead to the exact asymptotic behavior of the xc potential. 
Eq. (20) holds for all /(r) that vanish in the asymptotic 
limit and for which the third term of Eq. 0 . does not 
contribute to the asymptotics of the functional derivative 
i.e., under very general conditions. Further de¬ 
tails of the calculation, specifically regarding the question 
of the xc potential asymptotics in different spin-channels, 
are given in Appendix [B| 

The LMF is limited between 0 < /(r) < 1 and there¬ 
fore the asymptote is bound between | < 7,7 < 1. Con¬ 
sequently, the exact value 70 - = 1 can only be reached 
by setting /(r) = OVr, which corresponds to the trivial 
case of using EXX, “as is” or combined with a purely 
(semi-)local correlation functional. 

A different extreme case, /(r) = IVr, does not, as one 


could naively believe due to Eq. (|21|), lead to 7 ^ = 4- 


Here, we have to take the neglected first term of Eq. (18) 


into account again, and from this we see that 70 - actually 
vanishes. This is to be expected, since in this case the 
local hybrid reduces to a purely (semi-)local functional. 

Fig.[T] shows a numerical verification of the above ana¬ 
lytical considerations (see Appendix for numerical de¬ 
tails). It depicts the xc (KLI) potential corresponding 
to the local hybrid of Eq. in comparison with the 
asymptotic decay according to Eq. (20) and Eq. (21) for 


the carbon atom. An additonal curve indicates the exact 
— 1/r decay, which is clearly not reached. The xc poten¬ 
tial, instead of decaying with 7 ^ = 1 as one would intu¬ 
itively expect [53], approaches the predicition of Eq. (21) 


( 7 ^ = 0.716) quite rapidly. 

Fig. shows the xc energy density exc(r) for the 
same system in comparison to its correct asymptotic of 
— l/(2r). Clearly the xc energy density shows the cor¬ 
rect asymptotic, cf. Eq. We thus see that while the 
behavior of e^c can directly be controlled via the LMF 
in Eq. (§, the process of finding the local xc potential 
via functional differentiation leads to non-local evalua¬ 
tions of the LMF that decisively impact the potential’s 
asymptotics. 

A physically meaningful quantity closely related to the 
asymptotics of the xc potential is the highest occupied 
eigenvalue Sho- Table [T] shows —Sho compared to the ex¬ 
perimental IP for the carbon atom for different function¬ 
als, together with the corresponding value of ja of the xc 
potential from Eq. (21). 


The LSDA, as generally known, significantly underesti¬ 
mates the IP due to the wrong potential asymptotics and 
the inherent self-interaction-error. Using pure EXX with 
the correct asymptotic decay and no self-interaction-error 
leads to a much better prediction of the IP. When em¬ 
ploying a local hybrid with the LMF /t(r), the explicit 
dependence on the parameter c becomes evident: With 
growing c, the asymptotic value 70 - grows and the de¬ 
scription of the IP improves. Fig. sheds further light 
on the situation. It shows potentials of local hybrids 
which are all based on Eq. (§ but use different values of 
c. Growing values of c increase the amount of EXX and 
lead to an overall deeper potential. This explains that 
the eigenvalues become more negative. 

We thus see that while all of the local hybrids used here 
can (so far, see caveat in the next section) be thought 
of as being one-electron self-interaction-free, they show 
different potential asymptotics and their highest occu¬ 
pied eigenvalues predict the IP with significantly dif¬ 
ferent reliability. The relation between freedom from 
self-interaction, potential asymptotics, and physical in- 
terpretability of the highest occupied eigenvalue as the 
negative IP is therefore much less clear than intuitively 
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TABLE I: Comparison of the highest occupied Kohn-Sham 
eigenvalue —Eho to the experimental vertical IP |93j for the C 
(£ 4 ^) and F atoms (£ 44 ,) using different functionals. All values 
are in hartree. 


System 

functional 

7aho 

^ho 

exp. IP 

G 

LSDA 


- 

0.2249 

0.4138 


/t(r)(c = 

0 ) 

0.6098 

0.2740 



/t(r)(c = 

0.5) 

0.7162 

0.3067 



/t(r)(c = 

1 . 0 ) 

0.7678 

0.3302 



/t(r)(c = 

2.5) 

0.8441 

0.3688 



ft{r){c = 

5.0) 

0.8966 

0.3970 



/o(r) 


0.8309 

0.3530 



EXX 


1.0000 

0.4378 


F 

LSDA 


- 

0.3808 

0.6403 


/t(r)(c = 

0 ) 

0.5055 

0.3810 



/t(r)(c = 

0.5) 

0.6665 

0.4724 



ft{r){c = 

1 . 0 ) 

0.7390 

0.5269 



/t(r)(c = 

2.5) 

0.8365 

0.6060 



/t(r)(c = 

5.0) 

0.8971 

0.6570 



/o(r) 


0.7927 

0.5798 



EXX 


1.0000 

0.6779 




FIG. 3: Asymptotics of the xc potential Uxct(r) for the C 
atom along the ai-axis, computed with pure EXX and local 
hybrids using /t(r) from Eq. with parameters c=0.5 and 
c=5.0. Also displayed are the corresponding asymptotics for 
-y^iEXX) = 1, 7 t(c = 0.5) = 0.716 and 7 t(c = 5.0) = 0.897, 
and the complete potential in the inset. 


believed. This observation also calls for taking a closer 
look at the iso-orbital indicator (/(r) that is used in en¬ 
forcing freedom from self-interaction. This is the topic of 
the next section. 


IV. THE IMPLICATIONS OF ORBITAL NODAL 
PLANES 


As explained in the preceding sections, many local hy¬ 
brids and other functionals such as meta-GGAs rely on 
the function ^(r) tending to 1 to detect regions of space 
in which a single orbital shape dominates the density, 



-15 -10 -5 0 5 10 15 

X (a.u.) 


FIG. 4: The function g{r) = on the numerical grid for 

the G atom 


and then, e.g., correct for self-interaction in such regions. 
However, a first caveat that one has to take note of is that 
g{r) —1 holds for one-particle densities of ground-state 
character. This is a possibly far reaching restriction for 
the use of ^(r) because electron orbital densities typically 
have nodes, i.e., are not of ground-state character. 

As a particular example, consider a case where the 
HOMO is a p-orbital with an azimuthal quantum num¬ 
ber m and is expressed in spherical coordinates as 
‘7’ho(r) = i?(r, At the region where the den¬ 
sity is dominated by the HOMO, Tw(r) = = 

||V 7 /n(r)p —i|Vi?(r, 0)p. However, r(r) —>• 
5|V(/3ho(r)P = + rn?R^{r,6)/{rs\ii6Y). As 

a result, while for m = 0, T(r) identically equals tw(i’), 
for TO = ±1 this is no longer the case, but r(r) and tw 
approach each other asymptotically, as the m^-dependent 
term of r(r) decays to zero. 

One may counter-argue that this restriction is not 
so severe because in density functional construction the 
condition d(r) —>• 1 is mostly used to detect those re¬ 
gions of space in a finite system which are far from all 
nuclei, where the density decays nodelessly, and there 
T(r) —)■ TwCr) even in the example above. 

However, we show below that even in such regions the 
condition T(r) —?> Tw(r) can be violated. This leads to 
a second caveat about the reliability of the g{v) indica¬ 
tor. It is rooted in the existence of orbital densities that 
have nodal planes or nodal axes. Fig. illustrates the 
case. It shows g{r) evaluated in the (a;z)-plane for the 
carbon atom (see Appendix for numerical details, in¬ 
cluding grid setup). The density here was obtained using 
/t(r)(c = 0.5), but the density features relevant here are 
not sensitive to functional details. The important obser¬ 
vation is that g(r) approaches I in the asymptotic limit 
in every direction - except for in the vicinity of x = 0. 

The first step towards an understanding of this finding 
is to note that the z-axis is a nodal axis, being the inter¬ 
section of the nodal planes of the two HOMOs of carbon, 
which are degenerate and of p-orbital character. 

The consequences of the existence of nodal planes 
can be studied analytically. To this end we look at 




















a schematic density that is dominated by the HOMO 
<Pho(r), but also take the next lower lying orbital (/Sho-i (r) 
into account, i.e. n(r) ^ |i^ho(r)P + |y)ho-i(r)P- With 
this ansatz one finds 


and 


(V|^hoP + V|y^ho-iP)" 

Sdv^hop + |(/?ho-lP) 


r _ l(|v^hoP + |VvJho-iH- 


( 22 ) 


(23) 


These two terms combined and evaluated on or close to 
a nodal plane (denoted by —>■), where (/?ho —>■ 0, yield 

n.p. 


Tw 


IV ^ho—11 


r n.p. |V(/3hoP + |V(pho-l| 


< 1 . 


(24) 


Even though yiho vanishes on the nodal plane, its gradient 
still yields a finite value and keeps the function g{r) from 
approaching 1. 

Fig. 1^ shows that the deviation from 1 has a noticeable 
spatial extension of a few a.u. This raises the question 
of how well the use of the iso-orbital indicator ^(r) leads 
to freedom from self-interaction, as in some regions that 
so far have been considered as iso-orbital ones, e.g., all 
space far from the system’s center, self-interaction effects 
may not be eliminated fully when the indicator aberrates 
due to the presence of a nodal plane or axis. A different 
interpretation of Fig. would be to reconsider one’s ex¬ 
pectation of where iso-orbital regions are, or what they 
are. The traditional point of view has been that all space 
far from a finite system’s center is of iso-orbital nature. 
Fig. and Eq. ( [24| ) may be interpreted to show that this 
is not the case when the HOMO has a nodal plane/axis 
extending to infinity. From this perspective one might 
say that ^(r) does exactly what it is supposed to be do¬ 
ing, i.e., it indicates that the nodal plane region is not 
of iso-orbital character. Yet, also from this perspective 
Fig.i reveals a surprising finding, namely that even in¬ 
finitely far from a finite system’s center, the density may 
not be of iso-orbital character. 

The nodal plane observation also forces us to take a yet 
closer look at the central topic of this perspective, the 
potential asymptotics. Nodal planes can influence the 
asymptotics of a local hybrid’s xc potential in two ways. 
First, it has been argued that all orbital-dependent func¬ 
tionals show non-vanishing asymptotic constants in their 
xc potential along nodal planes of the highest occupied 
Kohn-Sham orbital that extend to infinity. This was first 
discussed in Refs. [94l [95] for the case of pure EXX, and 
the occurring shift was determined to be 


C(J - '^'X.cMfjCr (^^) 

with fixczo- = J ip*^{r)v^ccr{'r)(p^„{r)d^r and u^cia = 
/ (p*^{r)uxcia{'r)<fia{'r)d^r. The index denotes the 
highest lying Kohn-Sham orbital that does not show a 
vanishing spin-orbital density along the nodal plane of 
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X (a.u.) 

FIG. 5: The LMF /t(r) = ——’ evaluated with c = 
0.5, on the numerical grid for the C atom. 


the HOMO. Since Eq. (251 follows from the KLI (OFF) 
equation without referring to a specific functional, non¬ 
vanishing asymptotic constants on nodal planes of the 
HOMO are expected on rather general grounds. 

Second, the fact that g{r) —1 is not guaranteed on a 
nodal plane can also affect the potential. For the sake of 
clarity, we again discuss this effect for the specific exam¬ 
ple of local hybrids. When the LMF tends to zero on the 
nodal plane, i.e., /(r) —>• 0 and Eq. (7) is obeyed, then 

n.p. 

the non-vanishing constant of Eq. ([Ml is the only effect. 
An example for this case is the LMF /t(r) with a finite 
value of the parameter c. It is depicted in Fig. for the 
C atom density in the (a;z)-plane, and one sees that there 
are no asymptotic features. This is because the reduced 
density gradient in the denominator causes /t(r) to van¬ 
ish in the asymptotic limit, regardless of the occurrence 
of a nodal plane. The potential decays like —j^jr in all 
directions, but along the z-axis a nonvanishing constant 




n.p. 


(26) 


appears. This is shown in Fig. j^for /t(r)(c = 0.5) and 
the F atom. One can clearly see how i'xct(r) decays with 
7 ^ = 0.6650, but, instead of zero, approaches a constant 
of = 0.0244, in agreement with Eq. (251. 

A different situation occurs when /(r) ^ 0, i.e.. 


n.p. 

the behavior of the indicator function along a nodal 
plane/axis of the HOMO prevents the LMF from reach¬ 
ing its intended limit. This happens, e.g., for fo(r) or 
/t(r)(c = 0) and is depicted in Fig. again for the C 
atom. The occurrence of a nodal axis here very clearly 
affects the LMF. Since in this case Eq. 0 is violated 
in the direction of the z-axis, the previous derivations 
cannot be used to predict the potential’s asymptotic be¬ 
havior. However, we have numerically checked the xc 
potential’s behavior. On the nodal axis it neither tends 
to — l/r, nor to Co- — but rather tends to some other 
value. Thus, the nodal axis in this case has a very no¬ 
ticeable influence on the potential asymptotics, which is 
hard to predict a priori. 












FIG. 6: Asymptotics of the xc potential Wxct(r) for the F atom 
along the (projected) 2 -axis (denoted z*, see Appendix [A| for 
definition), computed using /t(r) with the parameter c=0.5. 



-15 -10 -5 0 5 10 15 


X (a.u.) 


FIG. 7: The LMF /o(r) = 1 - on the numerical grid 

for the G atom. 


V. CONCLUSIONS 

With local hybrid functionals serving as an explicit ex¬ 
ample we have argued that freedom from self-interaction 
in the sense of Eq. ([^ does not necessarily lead to the 
expected — 1/r decay of the local Kohn-Sham xc poten¬ 
tial. We have further argued that the ratio of the von 
Weizsacker kinetic energy density to the positive Kohn- 
Sham kinetic energy density, which is frequently used in 
functional construction for indicating iso-orbital regions 
and eliminating self-interaction effects in these, may not 
serve its intended purpose because it is very sensitive 
to excited state features such as orbital nodal planes 
that are present in Kohn-Sham orbitals that construct 
ground-state densities of many-electron systems. 

These findings have immediate and somewhat discom¬ 
forting consequences for the local hybrid approach. For a 
large class of functionals one has to accept that the cor¬ 
rect long-range xc potential simply cannot be obtained. 
This observation plays a role in explaining why it is very 


hard to construct a local hybrid that yields good bind¬ 
ing energetics and physically meaningful eigenvalues with 
the same functional form and set of parameters [53]. How¬ 
ever, the impaired relation between self-interaction and 
the xc potential’s asymptotics, and also the impact of 
nodal planes, stand in a context that is much larger than 
the local hybrid one. The iso-orbital indicator g{r) has 
been used in many functionals, not only local hybrids. 
Nodal planes are known to impact the exact exchange 
potential in surprising wavs |94[[95] . They have appeared 
here as a prominent feature in kinetic energy ratios, and 
we expect [53] that they are playing a much larger role 
in the exchange potential than has been realized so far. 
The observation that a one-electron self-interaction-free 
energy can go together with a potential that does not fall 
of like —1/r is not only a feature of local hybrids, but 
has also been reported for a “scaled down” version of the 
Perdew-Zunger self-interaction correction m- One may 
therefore wonder whether semi-local indicator functionals 
are in some sense incompatible with the fully non-local 
self-interaction correction that is achieved by EXX or 
full Perdew-Zunger-type correction approaches. It has 
also been pointed out recentlylHH] that Eq. 0 itself, 
which is the basis of the present definition of one-electron 
self-interaction, leads to questions when evaluated for or¬ 
bital densities, because Exc[n] is intended to be used with 
ground state densities, whereas orbital densities are ex¬ 
cited state densities. Further conceptual questions about 
Eq. Q relate to its inherent identification of orbitals with 
electrons and its unitary variance ElMlIinS]. The suc¬ 
cess of self-interaction corrections schemes that rely on 
Eq.(§ tells us that the equation is meaningful. However, 
the sum of the insights into its limitations that emerged 
over the years suggests that there is more to the question 
of self-interaction in density functional theory. 

While the above considerations point out areas that 
require further thought and work, one should also note 
that there have been developments in DFT that shine 
a bright light into the future. The concept of many- 
electron self-interaction [irm[T02] is not as straightfor¬ 
ward to use as Eq. ([^, but it avoids the conceptual 
questions that are associated with this equation. Range- 
separated hybrids yield the correct asymptotic potential 
and have proven to be a very successful concept, with¬ 
out being self-interaction-free |33J [331 ll03H116j . There 
have been successful functional constructions that can be 
seen as combinations of the local hybrid and the range- 
separation idea unidii]. Ensemble corrections m 
allow to extract information from functionals in an un¬ 
expected way, and can, e.g., further improve IP predic¬ 
tion. Finally, it has recently been shown [13(1] that a new 
type of a generalized gradient approximation can show 
features that were so far thought of as being associated 
only with exact exchange, such as step structures and 
surprising nodal plane features [96] . and understanding 
potentials in terms of xc charges has provided new in¬ 
sights |121H123] . Therefore, the battle against DFT’s 
old foe, the self-interaction error, and its surprisingly in- 
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dependent side-kick, the wrong potential fall-ofT, is far 
from being lost. 
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Appendix A: Numerical details 

We used the all-electron code DARSEC |124j for all calcu¬ 
lations presented in this perspective. This code exploits 
the rotational symmetry of diatomic molecules along the 
interatomic axis z, treating the azimuthal angle analyt¬ 
ically and thus effectively reducing the problem of solv¬ 
ing the Kohn-Sham equations to two dimensions. The 
equations are represented on a real-space grid of prolate- 
spheroidal coordinates. In such a coordinate system, the 
nuclear position(s) coincide with the focal point (s) of the 
grid located at z = ±i?/2, with R being the bond length 
of the diatomic molecule. This is the case also for cal¬ 
culations of single atoms: the position of the nucleus is 
not equivalent to the origin of the coordinate system, 
but is located at z = —Rf2, where R was set to be 
0.5 a.u. E.g., the C atom in our plots is centered at 
vc = {x,z) = (0, —0.25). The cc-axis is defined as per¬ 
pendicular to the z-axis, crossing the latter at z = 0, i.e., 
at a point being equidistant from the focal points of the 
grid (see Ref. |124j for details). 

In order to avoid numerical instabilities due to singu¬ 
larities in the Laplacian, the grid was chosen such that 
it does not include the actual z-axis, i.e. the interatomic 
axis. As a consequence, in this direction all quantities 
can only be plotted along a projected z*-axis, which takes 
into account all grid points that are closest to the actual 
z-axis. Since the discrepancy between the projected and 
the real z-axis decreases with increasing number of grid 
points, we made sure that the difference between z and 
z* is small by choosing sufficiently dense and large grids. 


Appendix B: The asymptotic decay of the 
exchange-correlation potential in detail 

In the following, we present considerations about the 
asymptotics of the xc potential in the spin channel that 
carries the global HOM O (g ho), as compared to the other 
spin channel (dho). Sec. [loused the condition that /(r) 
needs to vanish at a sufficient rate in the derivation of 


Eq. (20). In the present work, we investigated two pos¬ 
sibilities for the decay of the LMF. 

First, ft (r) for a finite value of the parameter c vanishes 
exponentially because /t(r) ~ ~ In 

this case, all individu al t erms in each functional deriva¬ 
tive, (see Eq. (17)), vanish exponentially in the 

asymptotic limit as well, except for the second term in 
Eq. (18). Eventually, this remaining term is responsi¬ 
ble for the reduced asymptotic decay of Eq. (20) due to 
the nonlocal evaluation of /t(r). Consequently, the xc 
potential in both spin channels decays with — 7 o-/r. 

However, a different picture occurs when evaluating 
/o(r) = 1— ■ This function decays much more slowly 

than /t(r) with finite c, as Fig (|^ shows for the car¬ 
bon atom. Consequently, not all terms in the functional 
derivative originating from /o(r) vanish individually and 
more detailed investigations are necessary. 

Defining K{y) = n(r)(ex(r) — e“(r)), the functional 
derivative in this case reads 




Nfj 




r - r' 


1 


+ V:^:,(r).V 
1 

2n3 (r) 

-I- Vn5 (r) • V 




5t{v) 


K{v) 


(HI) 


with - ^w(r') _ Sfojr) Tw(r) q., „ 

<5T(r) “ T2(r') “ <5rw(r) r(r) ' J^nereiore, Dorn 

s^(r) ^each the same absolute value in the 

asymptotic limit, but show opposite signs. Now, we have 



FIG. 8: Comparison of /o(r) and /t(r) with c = 0.5 for the C 
atom along the a;-axis. 
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to distinguish between the spin channels: If one looks 
aT”* a (^) channel that has the global 


at u 


HOMO, i.e. n(r) ~ crho(’')P> then one can see from 

Eq. (B1) that the fourth and fifth term are equivalent in 
the asymptotic limit except for the sign. Therefore, they 
cancel each other and, since the first and second term 
decay fast enough, only the third term remains, lead¬ 
ing to the limit of — 70 - 1 , 0 the other spin channel 
however, the fourth and hfth term do not cancel any¬ 
more, since the density is still dominated by oho('')j 
whereas the fourth term features g.j,„(r). Therefore, 
in the other spin channel yet another asymptotic limit is 
obtained, again strictly following from the evaluation of 


the functional derivative. 

This feature can be corrected by using a spin-polarized 
ansatz with an indicator function that is a spin-polarized 
LMF of the form 3 ^( 1 ") = with rwo-(r) and ^^(r) 

being the kinetic energy spin densities. In this case, a 
functional derivative that does not feature the total den¬ 
sity n(r) follows and therefore the aforementioned effect 
does not occur. However, since for the spin channel Uho 
all derivations made are valid independently of the form 
of the LMF and since this spin channel features the phys¬ 
ical meanigful quantity —Eho, it suffices for this work to 
consider the more simple LMFs instead of their spin- 
polarized counterparts. 
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